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Abstract 

We determine the quark mass ratio nic/ms on the lattice, using Wilson-type fermions. 
Configurations with Nf = 2 dynamical clover-improved fermions by the QCDSF collabora- 
tion are used, which were made available through the ILDG. In the valence sector we use 
a sophisticated, mass- independently 0(a)-improved Wilson- type action with small cut-off 
effects even in the charm mass region. After an extrapolation to the physical pion mass, 
to zero lattice spacing and to infinite box volume, we find rric/ms = 11.27(30)(26). 



1 Introduction 

Quark masses are among the fundamental parameters of the Standard Model of particle physics. 
As they cannot be measured directly, their determination involves a substantial amount of 
theory - for decades uncertainties have been hard to estimate and error-bars were large In 
recent years Lattice QCD has made enormous progress at pinning quark masses down with a 
few-percent accuracy; see e.g. [2] for a summary. For ratios of quark masses the situation is even 
better, since in this case no lattice-to-continuum matching factor (whose accurate determination 
represents one of the most demanding steps in such a computation) is needed. 

The charm-to-strange quark mass ratio mc/ms (which is scheme and scale independent) is 
of direct phenomeno logical relevance |3]. It has been determined by HP QCD P] and ETM |5]. 
Both collaborations use lattice formulations with small cut-off effects even in the charm quark 
mass region, albeit with isospin (or taste) symmetry breaking, i.e. the pions are non-degenerate, 
in spite of a single being used, an effect which disappears oc a? with a the lattice spacing. 
By contrast unimproved or 0(a)-improved Wilson fermions avoid such effects, at the price of 
having comparatively larger cut-off effects (see App. A of [2] for a discussion). 

In ^ we constructed a Brillouin-improved Wilson action which was claimed to show small 
cut-off effects without isospin breaking, thus allowing for a one-to-one identification between 
lattice and continuum flavor. The latter feature is important, as isospin breaking effects require 
a more involved analysis, rendering it less transparent. Here we test the smallness of the cut-off 
effects by calculating the ratio m^/ms in this formulation (with tree- level clover improvement 
and one step of link smearing) in the valence sector (for s and c). The lattices with 2 degenerate 
dynamical flavors (for u and d) are provided by the QCDSF collaboration. The remainder of 
this article describes how we calculate the ratio on each ensemble, and how we remove the 
lattice artefacts to find the physical value of mc/rris- We end with an illustration of how this 
ratio may be used, together with a precise m^ input, to yield a robust estimate of m^. 
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2 Strategy to compute mc/rris on each ensemble 

Our goal is to compute the quark mass ratio nic/ms with controlled systematics. We follow a 
two-step procedure. In the first step we tune, for each ensemble, the bare mass parameter n of 
our action (see [6]) to the physical strange or charm quark mass and evaluate rric/rns on that 
ensemble. In the second step we eliminate the lattice artefacts by means of a global fit. 
Our strategy to compute m^/ on a given ensemble can be summarized as follows: 

1. Tune and Kc at the same time such that M^J^Ml.-Ml^ ) and {2M^-Ml )/(M|,,-M|,^ ) 
take their physical values of 0.80138 and 12.402, respectively [1]. These numbers build on 
M^^ =0.6858(8) GeV for the quark-line connected state, which follows via {2M]^-M^y/^ 
with SU(2)-symmetric values of Mk,Mt^ from [2], or from a direct computation 

2. Determine for either tuned k the PCAC quark mass, and form the ratio r = m^^^'^ / m^^^'~^ . 
In this step m^^^^ is determined from the connected ss correlator, while m^*^"^^ follows 
from the cs correlator together with the strange mass determined before (see below). 

As a theoretical caveat let us remark that in general with Wilson-type fermions the sea quarks 
affect the renormalization properties of the valence flavors. For a bare PCAC quark mass [7j 

AWI _ -^^^PCACTi , ^i. !, N„™W 



1 + {bA-bp)amY + {bA-bp)aTT{M) + 0{a' 



where mj is the Wilson mass of flavor j, and M the quark mass matrix. The Zj with J G {A, P} 
are lattice-to-continuum matching factors, while bj = l+0{as), 6j = 0(ag) denote improvement 
coefficients. As we follow a tree- level improvement strategy (with csw = l) see [6]) the ratio 

1 + 0{a,a) + 0(0^)1 (2) 



^AWI ^PCAC 



™AWI ™PCAC 
k k ~ 

is found to carry two types of cut-off effects. As we shall see, the lack of knowledge which type 
would numerically dominate creates a major source of systematic error on the final result. 

Once r = m^^^^ / m^^^'^ is in hand for each ensemble, the final answer follows through three 
more steps (which, in practice, will be combined into a single global fit): 

3. Correct, for each ensemble, the value of r for the effect of the finite spatial volume L^. 

4. Extrapolate, for each /3, the result of step 3 to MP'^y'^ = 134.8 MeV [2] in the sea. 

5. Extrapolate the result of step 4 to the continuum, using an 0{asa) or O(a^) ansatz. 

To test how reliably the systematic uncertainties are assessed, we will repeat steps 3-5 for the 
control quantity M'i / {M})* — M"^ ), whose physical value is known. 



3 Analysis details and final result for rric/ms 

We now give details of how we determine the ratio rric/ms on each ensemble, and how we 
eliminate the lattice artefacts by means of a global fit. 

We use the Nf = 2 configurations by QCDSF [8HTT] made available through the ILDG [T2] . 
Since we measure dimensionless ratios, one might naively think that no scale determination is 
needed. However, in the extrapolation to the physical point a scale is required. We will use [Hj 

a[fm] = 0.076,0.072,0.060 at /3 = 5.25, 5.29, 5.40 (3) 
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Table 1: Details of the QCDSF Nf = 2 lattices made available through the ILDG, with aMj^ 
from fWI (in one case inferred from The values of M.^, L in the same block are based on 

the scales For comparison we add another block with information on Mt^, a, L from 



/3 


Ksea box size confs 


1/Ks 1/Kc O3 O4 


5.25 


0.13460 16^x32 2-500 
0.13575 24^x48 2-500 
0.13600 24^x48 2-500 


7.8310(18) 8.793(19) 13.48(22) 2.220(54) 
7.8504(12) 8.612(15) 12.87(11) 2.081(38) 
7.8520(11) 8.548(11) 12.89(12) 2.114(36) 


5.29 


0.13550 24^x48 2-400 
0.13590 24^x48 2-500 
0.13620 24^x48 2-500 
0.13632 24^x48 386 
0.13632 32^x64 2-500 


7.8601(13) 8.641(13) 12.69(16) 2.123(39) 
7.8632(15) 8.574(16) 12.68(16) 2.089(42) 
7.8635(09) 8.502(09) 12.85(14) 2.056(34) 
7.8615(17) 8.477(17) 13.41(25) 2.031(57) 
7.8648(08) 8.478(09) 12.57(08) 1.977(27) 


5.40 


0.13560 24^x48 2-500 
0.13610 24^x48 2-500 
0.13640 24^x48 2-500 
0.13640 32^x64 2-500 
0.13660 32^x64 2-500 


7.8823(10) 8.503(10) 13.01(15) 2.163(46) 
7.8859(10) 8.463(09) 12.55(17) 2.073(37) 
7.8842(11) 8.403(12) 12.66(19) 2.041(51) 
7.8864(08) 8.417(08) 12.25(11) 1.962(28) 
7.8862(08) 8.397(08) 12.42(11) 1.955(29) 



Table 2: Tuned kappas of the Brillouin operator and final 03^4 for each ensemble. Usually 500 
configs were downloaded; in most cases inversions were performed on more than one timeslice. 

for this purpose, but apart from the extrapolation this scale is not used. Given the resources 
available to us, we select the 13 ensembles marked with a bullet or circle in Tab. [1] for analysis. 
They cover a wide range of pion masses and box volumes (both in fm and in MttL units), so that 
a controlled extrapolation to the physical pion mass and infinite volume should be possible. 
On a given ensemble, for a few mass parameters 1/ks,1/kc, we determine the correlators 
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of four mesons (the pseudoscalar and vector with cs and ss flavor content). From these we 
form the observables d = Mj,,J{M^,^- Mj,,^), O2 = (2M|,g, - Mj,J/{M^,^-Mj,,^), and 
03 = {2mcs—Tnss) /rUss, where rriij denotes the average of the PCAC masses with flavor i and j, 
based on the improved symmetric derivative d(j){t) = [0(t— 2) — 8(f){t—l) + 80(t+l) — 0(^+2)] /12. 
For each observable a spline interpolation in 1/ks and 1/kc is constructed. The target value 
Oi = 0.80138 defines a line in the (I/k^, 1/kc) plane, and the same holds true for 02 = 12.402. 
The point where these two lines intersect defines the tuned set (I/k*, and the value of 

O3 at this point is the desired ratio r on that ensemble. The spacing in l/ns and 1/kc is chosen 
sufficiently narrow so that the uncertainty due to the interpolation is completely negligible. 
Since all of this is done inside a jackknife, the jitter of the crossing point is fully propagated 
into the statistical error of the tuned r, as listed in Tab.[2l For an illustration see |13j . 

Finally, we wish to correct for the systematic effect that the finite lattice spacing (a > 0), 
the larger-than-physical pion mass (M^r > M^^^^) and the finite spatial volume {L^ < 00) have 
on the measured iric/ms, by means of a global fit to our dataset. For each artefact, we invoke 
an extrapolation formula which is consistent with both theoretical expectations and the data. 
We shall consider several (reasonable) options for each effect, and treat the spread of these as 
the systematic error of the final result. The dominant cut-off effects may be proportional to 
Ostt (what theory suggests) or proportional to (what empirical evidence seems to prefer [B]). 
In the range of interest the dependence on m^^^ may be a quadratic or cubic function of Mj^. 
Finite volume effects may be proportional to Ki{M^L) / (M^^L) exp{—M^L) / (M^L)^^'^ times 
M^/(47rF,r)^ (as in the case of Mt,{L)/Mt, — 1), or just proportional to (as frequently used 
in the old lattice literature). By combining these forms we arrive at the 8 ansatze 

r^''^'''\a,M^,L) = b ■ \l + Sf^{a) + S^^g^^\M^) + e^''^ h^''\M^ , L)] (4) 



withi,j,fc G {1,2}, where /« = asa, f^'^^ = a\g^^'^ = Mlg^^'> = M'^, /iW = ^M^/L^ exp(-M^L), 
= l/L^. Note that this is the first time that we make use of the auxiliary scales (|3]); here 
we need to invoke them, since the coefficients c'^'^\ e^^^ are dimensionful quantities. 

The last point to be discussed is which ensembles are included in the fit. It turns out that 
the ensemble 5.29_0.13632_24^ x48 cannot be described by any of the ansatze; once we drop it 
all versions of @j yield consistently x^/d.o.f ~1. For r*^*'-''''^)(0, M^'^^^, 00) one finds 11.01(36), 
11.02(32), 11.24(33), 11.05(35), 11.39(26), 11.41(24), 11.60(25), 11.43(26), respectively, where 
the errors are purely statistical. To avoid underestimating the effect of the extrapolation, we 
need to include the spread among these 8 results as a source of systematic uncertainty. This 
yields rric/ms = 11. 27(30) (22), where the standard deviation of the distribution is used as the 
systematic error. Since some of the ensembles feature large pion masses and small volumes we 
use the cuts M^^ < 670, 900 MeV and/or L > 1.4, 1.7fm to check for any additional systematic 
uncertainties. This yields six independent options for our dataset. The center of this enlarged 
distribution (from the 8 ■ 6 = 48 analyses) is lower than the central value mentioned above, 
amounting to an additional systematic uncertainty of 0.14 which we add in quadrature to the 
previous one. This yields our final result 

mjm, = 11.27(30)(26) (5) 

in the continuum, at the physical mass point, and in infinite volume. 

To illustrate the procedure we present one of the 8 global fits - the (i, j, k) = (2, 1, 2) variety 
with O(a^), O(M^) and 0{1/L^) terms - in Fig.[TJ The data have been shifted by the effect of 
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Figure 1: One of the 8 global fits, namely r'-^'^'^^ with O(a^), 0{M^), 0{1/L^) correction terms, 
for the joint extrapolation to zero lattice spacing, physical pion mass, and infinite volume. 
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those terms which are not on display. For instance, in the continuum extrapolation panel 

plotdata(a) = data(a, M^, L) - fit(a, M^, L) + fit(a, Mf'^^ oo) (6) 

is shown as a function of a^, while in the pion mass extrapolation panel the last term reads 
"fit(0, M^, oo)" , and in the infinite volume extrapolation panel it is "fit(0, M^^^*^, L)" . Note 
that this affects only the presentation, not the final result 

To test whether our assessment of systematic uncertainties is true and fair, we apply the 
same analysis procedure to the observable O4 = M^/{M]j* — Ml)J. This gives 1.79(08)(12), 
which agrees perfectly with the physical value 1.7707 [Ij. This supports the view that our 
analysis procedure yields reliable estimates of the uncertainties in ([5]). 

4 From quark mass ratios to individual masses 

To give the reader an idea of what can be done with our result ([5]), we combine it with an aggre- 
gate value of rric to obtain an estimate of m^. For several precise results are available, which 
use either sum rule techniques or perturbative estimates of moments of current correlators. By 
contrast, computing rris directly on the lattice involves renormalization factors like the factor 
Za/Zp in ([1]) whose non-perturbative determination is technically quite demanding. Therefore, 
computing via ([5]) from rric offers the possibility to check the current best calculations of 
rUs (see [2] for an overview) without recurrence to Z-factors [3]. 

We now collect the current best estimates of the charm mass, which have a 1-2% error. The 
first result mc(3 GeV) =0.986(6) GeV [H] is based on the current correlator method on the lat- 
tice. The remaining ones are based on sum rules and experimental electron-positron annihilation 
cross section data, namely mc(3 GeV) =0.986(13) GeV [15], mdrnc) = 1.277(26) GeV [16], and 
mc(3 GeV) = 0.987(09) GeV |T7], respectively (for an examination of the uncertainties involved 
see in particular [IS])- Through standard 4-loop MS running, these results can be evolved to the 
common scale /i = 2 GeV, where they read mc{2 GeV) = 1.092(7), 1.092(14), 1.096(22), 1.093(10) 
GeV, respectively. A straight mean of the central values and of the systematic uncertainties 
yields the conservative average mc(MS, 2 GeV) = 1.093(13) GeV [HHIIj. 

Upon combining this input value with our result we arrive at the estimate 

m,(MS, 2 GeV) = 97.0(2.6) (2.5) MeV (7) 

which does not build on a renormalization factor. At this point we may continue by using the 
ratios rris/mud = 27.53(20)(08) and {md — mu)/{md + mu) = 0.381(05)(27) by the Budapest- 
Marseille- Wuppertal collaboration [I8l|19], where mud = {rnu+md)/2, to end up with 

m„d = 3.52(10)(09)MeV , m„ = 2.18(06)(11) MeV , = 4.87(14)(16) MeV . (8) 

Still, the precision reached is competitive in view of the global averages given in [2j. 

This concludes our illustration how the light quark masses can be obtained without recur- 
rence to renormalization factors, at the price of including perturbative information. 

5 Summary 

The goal of this note has been to calculate the ratio rric/ms, using our relativistic fermion action 
[6] in the valence sector, with a controlled extrapolation to zero lattice spacing, to physical sea 
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pion mass and infinite box volume. The only systematic effect which is not controlled is the 
quenching of the strange and/or charm quark, but this is the case in other state-of-the-art 
calculations [HIS], too, and there are good reasons to believe that the effect is negligible on the 
scale of the error in ([5]) (cf. the discussion in [2]). 

Our result ([5]) is consistent with the values mc/m<( = 11.85(16) by HPQCD j4j and 12.0(3) by 
ETM ^ (note that the spread among the entries in their Tab. 7 has not been propagated into 
their final error), with a slight tension at the level of 1.36a and lA7a, respectively. Though 
nominally less precise, our result serves as an important benchmark, since our formulation 
bears the unique feature that it is free of any lattice-induced isospin (or taste) breaking. The 
relatively mild slope in aga or as determined by our global fits and the small overall spread 
among the entries in the 0^ = mc/ms column of Tab. [2] support the view that the formulation [6] 
entails small cut-off effects up to the region of the physical charm quark mass. 

For illustration we combine our ratio (jS]) with an average of nic from [T3HT7] to obtain the 
value ([7]) of mg. While there are results on with a higher claimed precision (see e.g. [2] for 
a review), our computation is the only one which avoids both Z-factors and unphysical isospin 
breaking effects, and this renders the result particularly robust and reliable. 
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